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1. INTRODUCTION 

The study of Virtual Compton Scattering (VCS), e + p — ► e! + p' + 7, at CEBAF 
and MAMI (Audit et al. 1993) could provide valuable information on the structure of the 
nucleon, complementing the information obtained from elastic form factors, real Compton 
scattering, and deep inelastic scattering. In this paper we concentrate on the kinematic 
region where the final photon has low energy — i.e., below the threshold for 7r° production. 
As shown by Guichon et al. (1995) the low energy cross sections are parametrized by 10 
generalized polarizabilities (GP), functions of the virtual photon mass. Their evaluation 
requires the knowledge of the nucleon excited states. This sensibility to the nucleon 
spectrum can provide substantial insight into the non-perturbative aspects of the QCD 
Hamiltonian. 

Guichon et al. (1995) made an initial evaluation of the GP's so as to provide an 
order of magnitude estimate of these new quantities and to illustrate their variation as 
a function of the virtual photon mass. In that calculation we neglected all recoil effects, 
that is terms which go like the velocity of the nucleon. As a result of that approximation 
only 7 GP's were non-zero. 

In this paper we extend the calculations to include the recoil corrections which turn out 
to contribute only when the final photon is magnetic. We also study the relationship of the 
nucleon excited states to the GP's. We use the Non-Relativistic Quark Model (NRQM) 
to take advantage of its simplicity and the ready availability of its wave functions. Also 
the separation of the center of mass and internal motion greatly simplifies the calculation, 
making it analytically tractable. In principle, it is possible to use other wave functions 
but this can be prohibitively laborious and messy. The NRQM estimate should be a use- 
ful guide to the analysis of the soon available experimental data on the p(e, e'p)^ reaction. 

2. GENERAL FORMS OF GP'S IN TERMS OF CURRENT DENSITIES 

We first briefly outline the formalism for the definition and the calculation of the GP's. 
We refer to Guichon et al. (1995) for a detailed account of the problem as well as for the 
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notations and conventions. The hadronic tensor (see Figure 1) is defined by 

S(p x -p'-q') 



H^ B (q'm' s ,qm s ) = J dp x ^ 



(N(p')\J»(0)\X(p x )) 



E N (q') + q'-E x (p x ) 



(x( Px )\r(o)\N( P )) 



+(N(p')\r(o)\x( Px ))- 



-(X(p x )\nO)\N(p)) 



+ Hseagulh (1) 



E N (q) - q> - E x (p x ) 

where J M is the hadronic current, X the intermediate baryon excitations and H^ agull the 
contact term generally required by gauge invariance. The reduced multipoles are then 
defined according to 



+m > +L +M A l 1 I o \ 

m'm s M'M 



r p'L'M',pLM i i i 



(L'M',LM\S s)H p N L B ' pLM (q'm' s , qm s ). 



(2) 



with 



H 



p'L'M',pLM f 
JVB 



(</m>,) = (47T)" 1 / dqdq'V:(p'L'M',q')H^ B (q'm' s ,qm s )V p (pLM,q) } (3) 



where V v (pLM,q) are the charge (p = 0), magnetic (p = 1) and electric (p = 2) basis 
vectors defined by Guichon et al. (1995). 

When p, p' are equal to or 1 the GP's are defined by 



p(p'L', P L)S^ 



1 n (p'L',pL)S ( , x 



qiL' qL nb 



(4) 



J q'=0 



In the case of a virtual electric photon the analogous definition does not yield a GP 
with the usual photon limit as q — > 0. As explained by Guichon al. (1995), one must 
therefore introduce mixed GP's according to 



fjp'L'M',LM 
NB 



(q'm' s ,qm s ) = (An)- 1 [ dqdq'V*(p'L'M' , q') H N i B (q'm' s , qm s ) (^M +1 (g))\ 
J i 

(5) 



p(p'L',L)S^ 



-H 



(p'L',L)S 



q >L' q L+l A - l NB 



(6) 



q'=0 



where 3^m(<?) * s the vector spherical harmonic. The 10 independent GP's needed to 
describe the low energy regime are then the following, 



p(ll,00)l p(ll,02)l p(ll,ll)0 p(ll,ll)l p(ll,2)l 
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p(01,01)0 p(01,01)l p (01, 12)1 p(01,l)0 p(01,l)l 



(7) 



In the low energy regime the following excited states contribute: iV*(| , 1535), iV*^ ,1650) 



N 



* / 3 " 



■ 2 ' 



1520), JV*(| ,1700), A(§ , 1232), A*(± , 1620) and A*(| ,1700) 



In the NRQM the current density in Eqn.(|I]) takes the form 

(X(p x )\J°(0)\N(p)) = N oPx (p x -p) 



(X{p x )\J(0)\N(p)) = N 



Px+ P 

6m n 



Px(Px-p) + p x{p x -p) 



+ 



2m, 



■^x(p x -p) x (Px-P) 



Here p, P and S are overlap integrals of current operators. If one takes into account the 
factorisation of the cm. and internal baryon wave functions, they can be written in the 
following forms. 



Px(Px-p) = J dpdXe-^ {p X - pyX (j) x (p,\,J x ,M x ,T x ,T X )Q(j) N (p,X,m N ,T A 



(9) 



PxiPx - p) = l^f dpd xe^ x -py* 



>x(P> A ' J x, M x, T x , t x ) (i Va -i Va)Q <f> N (p, \ m iv> t at)> 



(10) 



X x (p x -p) = JdpdXe *Vl(Px pyX (f) x (p } \, J x ,M x ,T x ,T X )cr 3 Q(J) N (p,\,m N ,T N ), 

(11) 

where Q = (| + t|/2) is the charge operator of the third quark and cr 3 is the twice the 
spin operator of the third quark. 

With specific internal wave functions for the nucleon ((j) N ) and the intermediate exci- 
tations {(f) x) j one obtains explicit forms for the current density and hence the hadronic 
tensor. 
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3. CURRENT DENSITY AND HADRONIC TENSOR IN NRQM 



To calculate the integrals p(p x — P), P(Px ~ P)i an d ^(Px ~ P) f° r intermediate 
states X being N*{\~ , 1535), N*{\~ , 1650), iV*(§~, 1520), N*{\~ , 1700), A(| + , 1232), 
A*(±~, 1620) and A*(|~, 1700), we use the wave functions from Isgur and Karl (1978). 
The calculation is straightforward though fairly lengthy. We arrive at the following ex- 
pressions. 



— I- 



.V8n q 



9 a 



e-*'«*T N (-l)t +M *y/2J x + l 



1 



J 



x 



M x -m N m N -M x 



lMx—m, 



.(?)■ 



(12) 



'M x - mjv - 



(13) 



S iV(28)(9) 



( 1 1 7" 

1 2 J * 



(14) 



1 I Jx 



[Mx-fi fi -M x 



(15) 



(16) 



(17) 



(18) 
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Aq) = s , 

AT 

where x are spin wave functions, t n the isospin quantum number of the nucleon (i.e., 
±| for (£)), and e m is the mth component of the spherical basis vectors. Eqns.(|T^|-|T^) 
are for intermediate states with the same isospin quantum number as the proton. 

The main characteristics of these integrals are summarized below: 

a. ) The A(1232) and the 4 8 component of the iV*'s contribute only to E(p x — p) . 

b. ) The 2 8 component of the iV*'s, and the A*'s contribute to all p(p x —p), P(Px~P)^ 
and E(p x -p). 

c. ) For small x = \p x —p\ the behavior of the integrals are Ea(se), P{ x ) oc 0(1), and 
£jv*, a* oc 0(x). 

The leading term and the recoil term of the hadronic tensor in Eqn.(JT]) are separated 
in the following way. We work in the initial cm. system so that p = —q. In the 
direct term of Eqn.([|), p x = and q' = —p' . So, aside from the energy denominator, 
the direct amplitude factorizes into a product of a qr-dependent current density Jd,XN{o) 
and a q' dependent current density Jd,Nx(q'), where 



J d,XN 

(q) = (X(p x )\J(0)\N(p)) 

direct 



Jd,Nx(q') = (N( P ')\J(0)\X(p x )} dl rect = N 



q Px(l) + p x(q) + 7^^x(q) x q 



6m q 2m q 



q ' Px(q') + Px(q') - ^x(q') x 4 



Qm q ' 2m q 



The current density in the cross term has a more complicated q and q' dependence. 
Because p x = —q — q', p + p x = —q' — 2q and p + p' x = —q — 2q', it involves terms 
depending upon both q and q' . Let 

(X(p x )\J(0)\N(p)) cross = J c ,XN(q') + 5JxN(q,q'), 

and 

(N(p')\J(0)\X{p x )) cross = J c , NX (q) + +SJ NX (q,q'), 

then Jc,XN(q') , J c ,Nx(q), 5J XN (q,q') and 5J NX (q,q') are given by the following ex- 
pressions: 

JcXN(q') = N 



q Px (-q') + P x (-q') + ^x(-q') x (-<?') 



6m q 2m q 



6 



J c ,Nx(q) = n 



6m, 



-p* x (-q) + p* x (- q ) 



2m, 



-£x(-<z) x (-<?) 



5J X N(q, q') 



5J NX (q, q') 



Nn 



3m, 



-Pxi-q') 



Nn 



3m, 



-Px(-q)- 



From Eqn.(^|) we see that J° depends only on q or q', so that 5J° = 0. We define the 
leading term of the hadronic tensor H^ B {q 'm' s , qm s ) by neglecting the terms depending 
on both q and q' in the cross term, and the q' dependence of Ex{q + q') in the energy 
denominator of the cross amplitude. That is, 



^NB-LEADINciq m si q m s 



E 



Jd,Nx(q') Jd,XN(q) 
M-M x 



J£xM) J^Nx(q) 

E{q)-E x {q) 



1 seagull' 

(20) 



The recoil term of H^ B (q'm' s , qm s ) includes all the effects coming from terms depend- 
ing on both q, and q', and those arising from the expansion of the energy denominator of 
the cross term to order q'. The recoil contribution contains terms to order q' and higher 
and can be written as follows, 



H NB-RECOIL(q' m 's^q m s) = £ 

X^N 



JZxnW) Sr NX (q, q') + SJ XN (q, q') J» NX {q) 



[E(q)-E x (q)f 



q ■ q 

Ex{q), 



E{q)-E x {q) 



J'xAq') J u c ,Nx(q) 



(21) 



To get the GP we need only keep terms to order q' . So the recoil effects contribute only 
to GP's with a magnetic final photon, /i ^ 0. 

With Eqns. (|2Ti| , |2"T| ) for the leading and recoil terms, one can carry out the partial wave 
projection using the definitions of section 2. The partial wave projection poses no partic- 
ular difficulty except the need for careful book-keeping. We therefore skip the details of 
the partial wave decomposition and give the final expressions for the GP's in section 4. 



4. GENERALIZED POLARIZABILITIES IN THE NRQM 



Here we give the final results for the 10 GP's in the NRQM. We also study the 
properties of the GP's in relation to the parameters determining the nucleon structure in 
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the NRQM. Our aim is to find those properties of the nucleoli structure to which the GP's 
are most sensitive. This should help identify the most useful aspects of VCS in studying 
the nucleon structure. 

The analytic expressions for the 10 GP's can be written in the following form, where 
the leading contributions are the same as in Guichon et al. (1995). (Note, however, that 
the curve for p(° 1 ' 1 ) s ' in Guichon et al. has an error where the cross term was a factor of 
2 too large. We correct it here. The overall sign for pC 01 ' 1 ) 5 " is also changed to conform 
with the definition for the electric virtual photon case.) 

The GP's are plotted in Figure 2. 

4.1 Leading Contributions To The Generalized Polarizabilities 



p(oioi)S = 1 1 1 2 ( Z d Jx , Z c Jx \ (22) 

-2S + ll8a* e X= ^ A , X \M-M X + E(q) - E x (q) ) 1 J 



1 fo i o o ( iUx-1/2 / 7 2,S,J X V2,S,J X \ 

p(0112)l = 1 6 1 -q 2 /6a 2 y^ 2 { jJ ( f 'ad f >ac \ {t y* 

108V5m g a 2 x= j^ A , x 2I X \M-M X E(q) - E x (q) J 1 } 



7S 

'A 

para 2S + 1 27 mf \M - M A ' E(q) - E A (q) 



Pr ,5 = i5D^e^ (25) 
54 m q a z 

The mixed GP is the sum of two terms: 

p(01,l)S = pj01,l)S + p(01,l)S ^6) 

p(oi,i)5 = 1 1 ]_ e -<7 2 /6a 2 a 2 ( Z d Jx Z c' Jx \ m) 

F 2S + ll08m q a 2 x=n* a* X \ M — M x E(q) - E x (q) J V ; 



iii ( i Ua-l/2 / 7l,S,Jx 7l,S,J x 

p(01,l)5 = 1 1 1 -g 2 /6a 2 V- „2 l^ 1 ) / Z ad Z ac 



a 



25+ l36>/3nv* 2 X= ^ A » 24 V. M - M * E(q)-E x (q)J 

(28) 
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Angular functions Z ,Jx , Z L,s,Jx and in the above summation are given by: 
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Jx 






rvL,S,Jx 


ZL,S,J X 
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4.2 Recoil Contributions To Generalized Polarizabilities 



P. 



(1100)1 



recoil 



3V3 m q 



E 



2 7^ 
^X^llOO 



£x(g)[£(g)-£ x (g)] 2 L 



1 - 



£ x (g)[£(g)-£ x (g)]' 



3a 2 



(29) 



3(1102)1 

recoil 



1 1 

3^ m q 



-q 2 /6a 2 



E 



« 2 7 J 



102 



X=iV*,A 



3a 2 



(30) 



3(1111)0 

recoil 



1 a 2 q 2/6a 2 



3 m 2 



E 

X=N*,A* 



n 2 7 Jx 



E x {q)[E{q) - E x {q)Y 



(31) 



5(112)1 
recoil 



where 



1 1 



6V15 m 2 



E 

X=N*,A* 



„2 7 Jx 



£ x (g)[£(g)-£ x (g)] 2 



E x (q)[E(q) - E x (q)} 



3a 2 
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^1100 
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^1111 
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The 4 8 component of the excited state wavefunctions contribute nothing to the proton 
GP's, because of the isospin factor (1 — 2r). However, they do contribute in the neutron 
case, which we do not study here. The parameters m q = 350 MeV and a = 320 MeV 
are taken from Isgur and Karl (1978). The P-wave intermediate states are ordered as in 



9 



the following table, together with their representation mixing parameters ax for the 2 8 
representation. 



X 


N*{\ ,1535) 


N*{\ ,1650) 


N*(\ ,1520) 


N*{\ ,1700) 


A*(| ,1620) 


A*(§ ,1700) 


a x 


0.85 


-0.53 


0.99 


0.11 


1.0 


1.0 
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5. N* PROPERTIES AND THE GP'S 



The NRQM parameters are well determined by fitting the static properties of baryons 
(Isgur and Karl, 1978). Nonetheless, there are other phenomenological models than the 
NRQM, such as the bag models. Different models may not necessarily give exactly the 
same properties for the nucleon and its excitations. We do not intend to survey all 
the different nucleon models in this paper, but just to investigate those nucleon and N* 
properties which exert the most important influence on the behavior of the GP's in the 
NRQM. Two main factors are studied here: the mass spectrum of the nucleon and the 
size parameter a. The results are displayed in Figures 2, 3, and 4. 

The GP's have a strong dependence on the mass and energy spectrum of the excited 
states of the nucleon and the A. In Figure 2 we used the average masses of the N* and 
A* from the particle data tables (1994). However, these masses are all determined within 
a range and may be different from the predictions of the NRQM. We study the effects 
of the N* mass spectrum by comparing the GP's calculated with the lower and upper 
limits of the masses from the particle data tables (1994) and also with those predicted in 
NRQM of Isgur and Karl (1978). Figure 3 shows that some GP's are very sensitive to the 
iV* masses, particularly at small q. The effect on pt 01 - 1 ) is especially large as compared 
with the theoretical masses of Isgur and Karl (1978). The change in pt 11 - 00 ) 1 , pC 11 ^ 2 ) 1 , 
and pt 11 ' 2 ) 1 are quite drastic due to a more complicated factor of the mass and energy 
differences. 

The effects of the hadron size parameter a, are illustrated in Figure 4, where the plot 
is for a variation of ±5% in a from its normal value of 320 MeV. The main influences are 
again seen to be in the low q region. 

In conclusion, we have presented a calculation of the N* contribution to the generalized 
polarizabilities for virtual Compton scattering on the proton. The dependence of these 
GP's on the N* properties has been studied. We hope that there will soon be experimental 
data with which these estimates can be compared. 

This work was supported by the Australian Research Council. 
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FIGURE CAPTIONS 



Figure. 1. Direct, Cross and Seagull Terms of the hadronic tensor in the lowest order 
QED perturbation theory. 

Figure 2. The 10 GP's in the NRQM, The parameters m q = 350 MeV and a = 320 
MeV. Note that pt 11 ' 00 ) 1 , pt 11 - 02 ) 1 , and pt 11 - 2 ) 1 are all zero in the absence of the recoil 
correction. (Note that in Fig. 2c the superscript "para" refers to the paramagnetic con- 
tribution from the A(1232) and "dia" to the seagull contribution which has the opposite 
sign.) 

Figures 3. Effects of the masses of nucleon(A) excitations on GP's. Four groups of 
different masses for P-wave excitations are used. 

dotted line - using the lower limit of the masses from the particle data tables (1994): 
iV*(1520) iV*(1640) 7V*(1515) iV*(1650) A*(1615) A*(1670) 

dashed line - using the upper limit of the masses from the particle data tables (1994): 
iV*(1555) iV*(1680) iV*(1530) iV*(1750) A*(1675) A*(1770) 

dot-dashed line - using the theoretical masses given in Isgur and Karl (1978): iV*(1490) 
iV*(1655) iV*(1535) iV*(1745) A*(1685) A*(1685) 

solid line - same as Figs. 2, using the average masses from the particle data tables (1994): 
iV*(1535) iV*(1650) iV*(1520) iV*(1700) A*(1620) A*(1700) 

Figures 4. Effects of a ±5% variation (304 and 336 MeV) in a from its normal value of 
320 MeV. Other parameters remain the same as in Figs. 2. The dotted line is for a = 304 
MeV, the dashed line is for a = 336 MeV, the solid line a = 320 MeV, as in Figure 2. 
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Figure 1 
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Figure 4 



